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I Abstract 

In hep-th/0402219 we considered large N zero-coupling (i-dimensional U(N) gauge 
theories, with Nf matter fields in the fundamental representation on a compact spatial 
fNj ■ manifold S time, with Nf/N finite. This class of theories undergo a 3 rd order 

deconfinement transition. As a consequence it was proposed that the dual string theory 
has a 3 rd order phase-transition to a black hole at high temperature. 

In this paper we argue that the same conclusions are valid for such theories at any 
0^ ' finite order in perturbation theory in the 't Hooft coupling A and Nf /N . It is conjectured 

that this continues to hold at strong coupling and finite value of Nf/N < 1/A, which 
^ \ suggests that the supergravity approximation to the dual string theory has a 3 rd order 

thermal phase transition to a large black hole. 

o 



On 



1 Introduction 



The consideration of large N gauge theories, their confinement /deconfinement phase-transitions 
and the possible relation to the transition from a thermal ensemble to a large black hole at high- 
temperatures in the dual string theory is a topic which is being actively discussed [1]. An important 
insight [2]-[4] was made by considering the gauge theories on a compact space, as there is an addi- 
tional parameter RA which can be tuned to weak-coupling allowing for a perturbation expansion, 
where R is the size of the compact space, and A is the dynamical scale of the gauge theory. As 
the gauge theory is on a compact manifold, the Gauss law constraint induces interactions between 
gluons and the matter multiplets even at zero coupling. 

Sundborg [2] and Aharony, et al. [3, 4], have provided a general framework to discuss these 
issues in the context of a unitary matrix model. This strategy has been developed and extended 
in a variety of ways [5]-[8]. The particular system that concerns us in this paper is large N, U(N) 
gauge theory with Nf matter multiplets in the fundamental representation with Nf /N finite. 

We previously have shown [5] that large N, zero-coupling, U(N) gauge theory with Nf matter 
fields in the fundamental representation on a compact spatial manifold S ,rf_1 x S 1 , with (Nf/N) 
finite, undergoes a 3 rd order confinement /deconfinement phase transition at a temperature T c . The 
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low-temperature phase (phase Aq) has free energy of O(Nj), interpreted as a gas of color singlet 
mesons and glueballs, while the high-temperature phase (phase A\) has a free energy characterized 
by N 2 f(N f /N, T), interpreted as that of a gas of gluons and of fundamental and anti-fundamental 
matter states. 

An intriguing question is the identification of the string dual to the gauge theory and the 
associated deconfinement transition. Witten [1] and others [2]-[8] have associated the confine- 
ment /deconfinement transition with the Hawking-Page [9] transition of the dual theory. When 
(Nf = 0) the low-temperature phase has free energy of 0(1), while the high temperature phase has 
free energy of 0(N 2 ), signaling the deconfinement of color singlet glueballs into gluons. In the dual 
string theory this is associated to a transition from a thermal gas of supergravity modes to a large 
black hole [1]. At zero-coupling of the gauge theory on S ,rf_1 x S 1 , this transition is 1 st order at the 
Hagedorn temperature Tjj [2, 3]. For pure Yang-Mills theory on S 3 x S 1 , computed to 3-loops in 
perturbation this remains a first-order transition [3, 4]. At weak-coupling in the gauge theory, on 
expects the thermal string/black hole transition to occur for black holes with size of order of the 
string scale, while a comparison with the supergravity approximation requires strong-coupling in 
the gauge theory. 

When (Nf/N) is finite, at zero-coupling the transition is 3 rd order on the gauge side [5]. 
However, no example is yet known of a string theory with a 3 rd order transition to a black hole. 
The challenge is to understand if such a transition exists in the corresponding string dual. If 
the zero-coupling results of a 3 rd order transition survives extrapolation to strong-coupling, then 
one would have the (relatively) easier task of finding an appropriate dual in the supergravity 
approximation. This task itself is not without difficulties. For example, with d = 4 one might be 
tempted to consider a IIB model with N D3 branes and Nf D7 branes, with (Nf/N) finite in the 
large N limit. However, this simple picture is not possible as the number of allowed D7 branes is 
limited. One approach to this issue is to consider d=3, with N D2 branes and Nf D6 branes at 
finite temperature, including back-reaction [10]. Unfortunately the construction considered in ref. 
[10] only gave the high-temperature limit of the system, so there was no opportunity to observe any 
possible Hawking-Page transition. Another approach does allow for a back reaction in a d=4 string 
theory with Nf ^ [11], but again the putative phase-transition has not been discussed. Finally, 
if one has Nf probe D7 branes embedded in AdS$ x S 5 [12, 14], the Hawking-Page transition 
remains first-order even a strong-coupling but in the probe approximation [14, 6]. Perhaps this is 
not surprising as Nf/N — > in the large N limit in the probe approximation, so that this should 
be compatible with Nf = 0. Therefore the possibility of a 3 rd order Hawking-Page transition with 
Nf/N finite is unresolved. 

In this paper we continue the discussion began in ref. [5] where Nf/N is finite, but now with 
non-vanishing coupling constants. In Sec. 2 we consider the planar contribution at large N to S e g 
(the effective action) in perturbation theory, as a double expansion in the 't Hooft coupling A and 
Nf/N. This is presented in terms of a matrix model for U, given in terms of the zero-mode of the 
gauge field 

a = —[ Aq, (1.1) 
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U = e if3a = e ia/T 



(1.2) 



and 

u n = ^Tr(U n ). (1.3) 

We argue in Sec. 2 that at low-orders in the double perturbation theory, there is no solution 
to the saddle-point equations with u n = (n > 1). Since the eigenvalue distributions of U can be 
written as 

n=— oo 

where u n = u n and no = 1, thus there is no phase with p(9) = constant. 

In Sec. 3 we discuss in detail the phase structure as a function of temperature that follows 
from the S e s of Sec. 2 for the special case where (Nf/N) = 0(\ 2 ). In this example the saddle-point 
equations require solution of cubic equations, and already exhibit many of the main issues of this 
paper. There is a 3 rd order transition from phase Aq to A\ [5, 13] at a temperature T c below the 
putative Hagedorn temperature, a feature which appears to be a generic property of gauge theories 
with (Nf/N) finite. 

In Sec. 4 we adapt the strategy of ref. [8] to the issues of this paper. In eq'n. (4.1) we 
present the gauge invariant matrix model consistent with the double perturbation expansion in A 
and (Nf/N) of large N non-Abelian U(N) gauge theory. First we recover the zero-coupling results 
of ref. [5]. The saddle-point equations which follow from (4.1) are shown in general to exclude the 
phase described by p(9) = constant, to any finite order in the double perturbation series in A and 
(Nf/N). Thus in this case, the low-temperature phase is always in phase Aq, i.e., that described by 
p(9) 7^ constant, covering < 9 < 2ir. We argue, again in perturbation theory, that the transition 
to the single-gapped phase A\ is 3 rd order. 

Sec. 5 summarizes our results, and offers plausible conclusions regarding the string dual 
evaluated in the supergravity approximation. It is reasonable to expect that the supergravity limit 
of string theories dual to the class of gauge theories considered in this paper exhibit a 3 rd order 
transition from a thermal gas of supergravity modes to a large black hole. No such example is yet 
known. 

Issues not discussed in this paper are: 

1. small black-holes [8, 16], 

2. smoothing out the phase-transitions at finite iV [7, 8], and 

3. the logical possibility that at a temperature above T c , there exist phases A n such that p(9) is 
n-gapped [15]. 

Since even when (Nf/N) is 0(X 2 ), c.f. Sec. 3, the solution of the saddle-point equations require 
analysis of a cubic, the general analysis of the saddle-point equations can involve complicated, 
model-dependent algebraic equations. The principle conclusions of this paper do not depend on 
such considerations. 



3 



2 The Effective Action 



A. Perturbation theory 

Consider free U(N) gauge theories on S d ~ 1 x time, with Nf matter multiplets in the funda- 
mental representation of the gauge group with Nf/N finite as N — > oo. The partition function on 
this compact space (therefore subject to Gauss' law constraints) is a matrix model 



2Z(x) 



[dU] exp \J2 ~ \z B (x n ) + (-l) n+1 Z F (x n j\ ■ \tr{U n ) + tr(ET 



(2.1) 



+ J2~ \z B (x n ) + (-l) n+1 z F (x n )]tr(U n )tr(U +n ) 



n=l 



where zb{x) and zf(x) are the single-particle partition functions for the bosons and fermions in 
the adjoint representation, while Zb(x) and Zp(x) are the single-particle partition functions for 
the bosons and fermions in the fundamental representation, where 

x = e -P = e~ l ' T (2.2) 

for S*^ 1 with radius R=l. [R may be restored from dimensional considerations as needed.] Recall 
that 

u n = 1 Tr(U n ) (2.3) 

Now consider the same theory with non- vanishing gauge couplings g 2 , with A = g 2 N the 
't Hooft coupling. It is useful, in this section without a specific theory in mind to keep any scalar 
self-couplings and Yukawa couplings at zero. The partition function can be written as 

(2.4) 



Z(x) = J du n du n exp [-%(«„) 

The planar contribution at large N to S e ff in perturbation theory, in a double expansion in A and 
Nf/N, consistent with gauge invariance, is 

N~ 2 S ef f(u n ) = S e ff(u n ) 

= {]T G n {x)\u n \ 2 

n 

+ A ^ Gm,n{ X )( U mU n U- n -m + CC.) 
m,n 

+ A ^ ' G m)Jl) p(x) (u m u n UpU— n — m —p + c.c.) 

m,n,p 

+ ■■■} 

+ c.c] (2.5) 

^ ' n m,n 

+ A ^ ' Fn,m,p{x) [Un'U'm'U'p'U'—n—m—p C.C.] 
n,m,p 

+ A ^ ' Fn,m,q,p(%) [Un'Um'U'p'U'q u —n—m—p—q ~\~ C.C.] 



n,m,q,p 



+ 
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Consider the Feynman diagrams of the theory, where each insertion of a fundamental loop into 
an adjoint line of a diagram gives a factor of X(Nf/N), while a radiative correction by an adjoint 
line to a fundamental propagator gives a factor of A. [Since we are keeping scalar self-couplings 
and Yukawa couplings at zero, these are the only corrections to be considered in this section.] In 
(2.5) the leading powers of A and of (Nf/N) in the expansion of S e ff are exhibited explicitly. The 
various coefficients in (2.5) are a double power series in A and (Nf/N), beginning at A = and 
(N f /N) = 0. 

Near the saddle point of (2.5) one may integrate out the u n (n > 2), since at the critical 
point x c in phase Aq, \u n (x c )\ « |ui(x c )| for n > 2. Therefore the effective action for u\ near the 
saddle-point for u\ is 



+ [ -rf J {(ui+«i) [-C + D\ Ul \ 2 + E\ Ul \ 4 + 



N- 2 S eff { Ul ) = \a\ Ui \ 2 + B\ Ui \ 4 + ... 

+ (uj + ufj [F + G\ Ul \ 2 + H\U!\ 4 + ...] 
K + L\ Ul \ 2 + M\ Ul \ 4 + ...]} 



+ («? 
+ ... 



3 + uf 



(2.6) 



The leading order contributions to the coefficients in (2.6) are 

A = 0(1) ; B = 0{\ 2 ) ; 

C = 0(1) and positive at leading order ; 

D = 0(A 2 ) ; E = 0(A 4 ) ; 

F = 0(A) ; G = 0(A 3 ) ;H = 0(A 5 ) ; 

K = 0(A 2 ) ; etc. 



(2.7) 



Since each insertion of a fundamental loop in an adjoint line is accompanied by a factor of 
\(Nf/N), higher-orders in (Nf/N) result in higher powers of A compared to those displayed in (2.6), 
while the 0(Nf/N) term is that of a single fundamental loop with adjoint radiative corrections. 
Consider the schematic structure of the saddle-point condition which follows from (2.6) 

N -2 d&eff 

dui 

+ u l {A + ...) + o(\ (^)'j =0 

where the curly brackets are a power-series in u\ and u\, with coefficients a series in A. [We have 
only exhibited the leading terms explicitly in the curly brackets.] Recall that C is a power series 
in A, with its leading term 0(1) real and positive, while A has leading-term 0(1), F of 0(A), and 
B of 0(A 2 ). Thus, u\ = is not a solution to (2.8) to the order we are working. Further, since 
corrections to the constant term in (2.8) are a power series in A and (XNf/N), the general structure 
of the saddle-point equation in perturbation theory, emphasizing the constant terms, is 
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N _ 2 dS e ff 



dui 

n+l 

v x / 

n=l,m=0 



-_(fi) c+ f t *«(^ 

n=l,m=0 

+ ui{2£ «? + ...} 

+ + = 0. (2.9) 

The curly brackets represent a perturbative expansion in u\ with the leading term exhibited ex- 
plicitly, where the C n are independent of u\ and iii, but a power-series in A. Since C = 0(1), and 
positive at leading order, cf. (2.7), u\ = is not a solution to (2.9) in any finite order of the double 
perturbation theory, as the powers of A and (Nf/N) do not match. A solution u\ = to (2.9) 
would require 

oo / AT \ n 

E A " +m f-jvf) C„, m = C (2-10) 

n=l,m=0 V iV / 

which is not possible in any finite order of the double expansion, since C at leading order is 
independent of A and (Nf/N). If (2.10) has a solution, it must be a non-perturbative or all-orders 
result. As a consequence, u\ = does not solve the saddle-point equations in perturbation theory 
if (Nf/N) finite is considered as an expansion parameter of the theory. 

B. Phases 

The eigenvalue distribution of the matrix U can be written as 

-. oo 

n=— oo 

where U- n = u n and no = 1. The u n , with n > 2 can be expressed in terms of u\ near the saddle- 
point by means of the saddle-point equations that follow from (2.5), which is how (2.6) is obtained. 
Schematically, (not being explicit about coefficients and not distinguishing between u n and u- n ), 
one has 

n 2 ~ (-^)+x(-^)u 1 + Xul + ... 



N J \ N . 

- ; " ( + A//;» 2 + A (^jfj U! + ... 

+ X 2 ul + ... , (2.12) 

which displays the leading orders in A and (Nf/N) in (2.12), and similarly with its generalization 
to u n . Eqn. (2.12) reiterates that u n ^ are the only solutions to the saddle-point equations within 
the context of perturbation theory. 

Thus, for (Nf/N) finite in the large N limit, (2.11)-(2.12) and the saddle-point equations 
require, in the double perturbation expansion that 

p(6) / constant . (2.13) 
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Therefore, the allowed phases (in the notation of ref. [5]) are 

(i) (Aq) : p(0) completely covers the circle, but is not constant. (2-14) 
(ii) (A 1 ): p(9) = 

for a single finite gap in the distribution. The phase (Aq) is that described by the saddle-point, 
while (Ai) can be described by the methods of refs. [2, 3, 5, 8, 13]. The criterion that p(9) > 
implies from (2.12)-(2.13) that 

N I -\ + °(jf) + °^ + ° ( A jf) + • • • ( 2 - 15 ) 

in the double perturbation theory. The (Aq, A\) phase boundary occurs when p(0) = at a single 
point, which then implies that (2.15) is an equality. The transition to phase A\ is third-order, as 
can be deduced from an analysis analogous to that of ref. [5]. [It is a logical possibility [15] that 
there could be additional phases A n , with n finite gaps (n > 2), whose phase boundary occurs at 
temperatures above that of the (Aq,A\) transition, but this is outside the scope of our discussion. 
Further, this does not occur in theories with A = 0, (Nf/N) / [5].] 



3 Phase Structure of an Example 

In this section we discuss in detail the phase structure for the case when (Nf/N) = 0(\ 2 ). The 
purpose of the example is to illustrate a detailed discussion of the saddle-point equations, and 
evolution of the phase-structure with temperature. Although there are only two phases in this case, 
the details of the evolution with temperature depend on the relative magnitudes of parameters. 
When (Nf/N) = 0(X 2 ), insertions of fundamental loops into the diagrams of pure Yang-Mills 
theory can be neglected to the order in A 2 we are considering. Further, this example is interesting, 
as one can additively combine the contribution of fundamental matter to S e f f with the results of 
ref. [4], where pure Yang-Mills theory in d=4 was analyzed to three-loop order. 

A. General considerations 

Consider (2.6)-(2.7) correct to 0(X 2 ), with (N f /N) = 0(X 2 ) 
(2.6) have coefficients of that order so that to leading order 

N- 2 S e ff( Ul ) = ml\ Ul \ 2 + b\ Ul \ A 
-d(ui + ui) 

where 

m 2 = K(T H - T) , 

with 

K = 0(1), (3.2) 

b = 0(X 2 ), d > and 0(N f /N), (Nf/N) = 0(X 2 ), and where T H is defined as the temperature 
where m 2 = 0. Separate u\ = u\ + iu\ into real and imaginary parts, so that the saddle-point 
equations that follow from (3.1) are 

m\u\ + 26 u\ \u\ \ 2 — d = 

m\u\ + 2b u\ \u\\ 2 = (3-3) 



Only the first three terms of 



(3.1) 
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which requires u\ = 0; consistent with d = (Nf/N)C real. To the order we are considering, C is 
computed at A = 0, so that it has the same value as in ref. [5]. The restriction (2.15) now becomes 
Kl < 1/2 + 0(A) in phase Aq. 



In order to study phase Aq, we need the saddle-point solution to 
N~ 2 S ef f(ui) = m\ \ui\ 2 + b \ui\ 4 - 2 dux 

which requires 

26 u\ + m\ ui - d = . 
Solutions to (3.5) are characterized by the discriminant 



D 



27 



d_^ 2 
2b, 



(3.4) 
(3.5) 

(3.6) 



When D > there are 3 real roots from (3.5), while if D < there is 1 real root, and 2 complex 
roots which are not relevant, as u\ is real. If 



(3.7) 



d 


2/3 


o 




>> 


2b 


36 



then D < 0. 
If 



d 


2/3 


2 

mi 




<< 


26 


3^ 



(3.8) 



the sign of D depends on the sign of (mf/6), in which case the solution for u\ is approximately 
that of ref. [3, 4] except instead of u\ = 0, now 

Ml ~ (d/mf) > for m\ > . 

Since (d/2b) ~ 0(1) in this section, (3.8) always occurs for T sufficiently below Th- For example, 
m\ ~ 0(A) is sufficient for this to occur. 

In the next subsection we concentrate on a toy model, with 6 < 0, recalling that 6 < for 
pure Yang-Mills theory in d=4 [3]. The analysis for 6 > is qualitatively similar, which we do not 
present here. 

B. Toy model 

In the general / for n > 1, which leads to the constraint (2.15). In this section 

(N f /N) ~ 0(A 2 ), so that 

|«i| < 1/2 + 0(A) 

in the saddle-point phase Aq. One can simplify the analysis without sacrificing any of the qualitative 
features of the physics, so that for the remainder of this section we consider (3.4), but with u n = 
for n > 2. Then p{9) = 1/2-7T (1 + 2u\ cos 6), which implies that 



M < 1/2 



(3.9) 
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with the (Aq, A\) phase boundary occurring when \u\\ increases to \u\\ = 1/2 in the Aq phase, 
with the transition to A\ of 3 rd order. 

C. b < 

We track the evolution of the system for increasing T, beginning with T < Th- 
(i) ml > ; D > 

If D > 0, this requires 



|6|< 



2 (m{) 



2\3 



27 d 2 



(3.10) 



and an approximate solution to (3.5) is 



u\ ~ d/m\ 

2_ d_ 
27 H 



< 



1/3 



(3.11) 



[The two other real solutions to (3.5) are that of maxima of S e ff(ui), and thus not relevant.] For 
sufficiently large m\ \m\ ~ 0(A) is sufficient], \u\\ < 1/2 so that the system is in phase Aq, and 



The free-energy is 
so that 



N- 2 S eff ( Ul ) ~ -d 2 /m\ . 
FIT = S eff ( Ul ) 



F/T = -N 2 {d 2 /m\) 

t2\ 



0(Nf)/(T H - T) 



(3.12) 
(3.13) 

(3.14) 



since d ~ O(NffN). Thus (3.14) has the same form as (1.1) of ref. [5], and the interpretation 
of phase Aq as that of a gas of gluons and color singlet mesons. Therefore for sufficiently low 
temperatures one is in phase Aq. 

The saddle-point value of u\ increases as T increases. Consider 

(ii) D > 0, near D = 0; m\ > 

Using (3.6) we have 



1/3 



+ 5 m 2 , 



(3.15) 



and 



Ul 



d 

d 

2ml 



1/3 



1 - o 



5m 2 



1 - O 



5m? 
m\ 



(3.16) 



d S 

We verify that d ^ ff («i) > 0, so that this minimum is stable. 
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If u\ < 1/2, the system is still in phase Aq at this temperature, while if u± > 1/2 the system 
has already made a transition to phase A\ at a lower-temperature. 



(iii) D = 0; ra\ > 



d 

W\ 

d 



1/3 



2(m?) 



However, 



and 



d 2 s, 



eff 



(«i) = o 



a 3 5, 



e// 



(«i) = -|6| ui <0 



(3.17) 



(3.18) 



(3.19) 



so that (3.17) is not a stable minimum. The system evolves to \u\\ = 1/2 at which value it makes 
a transition to phase A\, if it has not already made this transition at a lower temperature. [See 
the discussion for (ii) above.] 

(iv) m\ < ; D < 

The solution to (3.5) is a single maximum at u\ < 0, and not a minimum. There is no stable 
saddle-point for |ui| < 1/2, so that at T = Th (i.e., m\ = 0), the system is already in phase A\. 

In summary, for (Nf/N) ~ 0(\ 2 ) and T sufficiently small, the system is in phase Aq. At 
some T c < Th there is a third-order phase transition to the high temperature phase A\. The free 
energy in phase A§ behaves as 



and in phase A\ 



F/T ~ Njfi(T) 



FIT ~ N 2 f 2 , T 



< T <T C 



T > T r 



(3.20) 



(3.21) 



The low-temperature phase is interpreted as a gas of (color singlet) mesons and glueballs, which 
makes a third-order confinement-deconfinement transition to a gas of gluons and fundamental and 
anti-fundamental matter states at high temperatures. 



4 Generalization 

We use the strategy of ref. [8] to discuss the issues of this paper suitable for all orders of perturbation 
theory. In doing so, we modify their set-up to suit our problem. Recall that u n = 1/7V Tr U n , so 
that (2.5) is a special case of 
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S(U,U + ) = J2 



a n (x) 



n=l 



n 



(Tr U n )(Tr U +n ) 



n=l 



b n (x) 



n 



(Tr U n ) + c.c. 



+ N 2 J2^k,-Ax)Y- k (U)Y- k ,(U + ) 



k,k' 



where 



^(3- krk ,(x)Y k (U)Y hl (U + ) 

k,k' 



Y k (U) = H(Tr(U/Nyp 

3=0 



[(Tr U n ) + (Tr U +n )] (4.1) 



?1=1 



(4.2) 



and where k and k' are arbitrary vectors of non-negative entries. Reality of the action and pertur- 
bation theory require [8] 

a k,k' = a t',k = a k',k ( 4 - 3 ) 
and similarly for f3 k k , . It is useful to define 



\k\ = J2 kj 

j 



(4.4) 



There is considerable redundancy in (4.1) which we remove in part by setting 

a k k = for \k\ = 1 



a 



k,0 



for \k\ = 1 



(4.5) 



0o,o = 



Contact can be made with the formulation of (2.5) by noting that the leading behavior in 
perturbation theory of the coefficients in (3.1) are 



a n (x) 


= O(l) 


b n (x) 


= 0(1) 


a k,k> 


~ (A)l*l+l*'l-2 


Pk,k' 





(4.6) 

Beyond the leading order of (4.6), the coefficients are a power series in A and (Nf/N), as well as 
any other coupling constants present in the problem. 



We take advantage of the Gaussian trick extensively exploited in ref. [8] to write 
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exp \J2 — tr(U n )tr(U +n ) 



,n=l 



P b r 

+ N f 2 — ^(^ n ) + ^(^ +n ) 



n=l 



n 



- g)7n(^)^e XP -^t(^ 

\ / n=l n=l 



HpJVj;^ [c n tr([/") + c n tr(U +n ) 



9n9n 



n=l 



where 



a n 9n + Jj- O n 



(4.7) 



(4.8) 



Eq'ns. (4.7)-(4.8) are closely related to (3.1) and (3.5) of ref. [5], except that now we allow a n and 
b n to be power series in (Nf/N) and A, with the leading order of a n and b n of 0(1), given by the 
zero-coupling expressions of (2.1). 

Define the free-energy F(g n ,g n ) by 



exp 



which implies that 



N 2 F(g n ,g n )] 

/( oo 
[dU] exp I N ]T ^ [tr{U n ] + tr(U +n ) 
K n=l 



1 °° 1 

F{_9n-> 9n) = T / , C n C n 
4 , n 

n=l 



(4.9) 



(4.10) 



in phase Aq. 

A. The Free-theory Revisited 

It is instructive to consider further (4.1), but first with y = fa y = 0. Using the Gaussian 
identity [8] 



exp 



71=1 ' 



n=l 



exp <J -N 2 X] ( ~ ) MnAtn + ^ 2 ( ) [^n#n + Hn9n] 



(4.11) 



The partition function is 



f N 4 V f p 
% = J J II dg n dg n dfi n dn n exp(N 2 S eff ) 



(4.12) 
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where here 



5 e// = - XI ( ~) + i X ( ~) ^ n 9n + Vn 9n] 

n=l ^ ' n=l ^ ' 



+ F(g k ,g k ) . 



The saddle-point equations are 



dS, 



eff _ 



dgk 
dSeff 



¥k + 



N f . . 
a k gk+ °k 



= 0; 



= - (j^J W + igk] = o 



(4.13) 



(4.14) 



(4.15) 



and two analogous equations for variations with respect to gk and yu.fc . Since ctfc(x) / 0, one has the 
solution 

(N f /N)b k 



9k = 9k 



(1 - a k ) 

since a k (x) and 6^(x) are real. The density of eigenvalues in the phase Aq is 



p{6) = — \ 1 + [fit exp(ifc6») + # fc exp(-z fc<9) 



fc=i 



so that from (4.8) and (4.16), one has 



Pk = 9k = Cfc/2 . 



(4.16) 



(4.17) 



(4.18) 



Thus, (4.16) is identical to eq'ns. (3.9) of ref. [5], where there is a 3 rd order GWW phase-transition 
[13] to phase Ai, interpreted as a confinement-deconfinement transition. Note that since (Nf/N) / 
0, g k = is not a possible solution, so that the phase p(9) = constant does not occur. Further 
discussion of the consequences of (4.9) and computation of p(9) for phase A\ may be found in ref. 
[5]. 

Thus, we have verified that the above formalism, with a k y = fl kk i = 0> reproduces ref. [5]. 
We now return to the general case described by eq'ns. (4.1)-(4.6). 

A. The General Case 

Equation (4.13) generalizes to 



Pnpr, 



S eff = ~J2 
n=l 

+ 1 X) ( ~) \i^9n + Un9n] ~ F(g k ,(j k ) 

71=1 

+ Xa^H) |fcWfc,| n(A)n,(/z) 

k,k> 

+ (^)E^(-o |S|+l ^'n(A)n'(M)EH) B i^+/in] 



(4.19) 



k,k' 



n=l 
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where 



and 



n^) = Ete)** 



(4.20) 



The saddle-point equations at large N that follow from (4.19) are 



as. 



dgk 



ff _ ( a k 



a l, Hi, + [ ^ ) b,, 



= 



(4.21) 



dS, 



eff _ 



J 



Hj + i gj 



+ E«M'H) |fc|+|fe,| ^n(M)n'(M) 

k,k' ^ 



+ 



+ 



^ ' ^ E/%,if(-0 |S|+|S ' 1 ^ EH) B (^n + An) 

k,k' 3 n=1 



_1 

iV 



N 



(4.22) 



k,k' 



and two analogous equations for dS e ff/dgj and dS e ff/djJLj. Since ctj / 0, (4.21) implies that 



= i 



"j ' ! - + 1 ^ ) 



(4.23) 



and similarly /ij = i[cij gj + (Nf /N)bj]. One then substitutes (4.23) into (4.22) to obtain an equation 
for gj in terms of a,j, bj, a^k' an d Pkk'i anc ^ analogously for gj. The solution to (4.22)-(4.23) for 
g~j and its analogue for gj determines the distribution p(8) in the ungapped phase, via (4.17) which 
is still valid. 

We now ask whether p(6) = constant is possible when (Nf /N) ^ 0. If so, this would require 
that g~j = for all j in (4.22)-(4.23), and analogously for gj. If k'j ^ 0, gj = gj = would require 



^ 6 



,\k\,\k'\ v 3 ' v 



Ifcl.lfc'l 



E(- 1 ) ri+lfc - 



* E ^ 

IfcMfc'l 



jvyxl^l+l*'! 



H) J+1 kfcW^w. 



(4.24) 



Note that cij and bj are both non-zero and of 0(1) in (Nf/N), while the right-hand side of (4.24) is 
a power series in (Nf/N), beginning with leading terms of 0(Nf/N) [see (4.5)]. Thus to any finite 
order in (Nf/N) (4.24) can not be satisfied by matching powers of (Nf/N). One concludes that 
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p(9) = constant is excluded in any finite order in the double perturbation theory. At low temper- 
atures the system is in phase Aq, and is not described by p(9) = constant at any temperature. 

The transition from phase Aq to A± is driven by the third-order phase transition in F(gk,gk), 
c.f. (4.10) and (4.13), which occurs when p(ir) = 0. This occurs when 



[Equality occurs if one may neglect for k > 2. See Sec. 3 for example.] Note that a,j(x) and bj(x) 
are continuous functions of x as one passes through the (Aq, A\) phase transition. From (4.23) this 
implies that fij and g~j(x) [and similarly Pj(x) and gj(x)] are smoothly related to each other at the 
(Aq, A\) boundary. Equation (4.19) can be written as 



where F(gj t ,gj t ) has a third-order transition at the (Aq,A\) boundary. Therefore it is plausible 
to assume that f(gk,gk) is smooth, with continuous derivatives at the phase-boundary. In the 
absence of an "accidental" third-order transition in f(gk,gk) which cancels that of F(gk,gk), or a 
lower-order discontinuity which seems unlikely given the smooth structure of f(gk,9k), S e ff has a 
third-order phase transition at the (Aq,Ai) boundary. 

5 Conclusions 

The phases of large N gauge theory at finite temperature on S d ~ 1 x time are described by p(0), 
the density of the eigenvalues of U in the unitary matrix model. When (Nf/N) is finite in the 
large N limit, p(0) = constant is not a solution of the saddle-point equations at any finite order 
in the double perturbation theory. The saddle-point equations describe a phase Aq, where p(9) ^ 
constant, but covers the interval < 9 < 2tt. There is a transition to a second phase A\, where 
p(9) is single gapped; the transition taking place at a critical temperature T c , which is lower than a 
putative Hagedorn temperature. The above conclusions are correct to any finite order in a double 
perturbation series in the 't Hooft coupling A and in (Nf/N). In order for the low-temperature 
phase not to be Aq, one presumably would require a phase-transition as one increased (Nf/N) at 
fixed A, or increased A at fixed (Nf/N), keeping the temperature fixed. We speculate that this is 
not the case. 

It is a reasonable assumption that generically f(gk,9k), c.f. (4.26), is sufficiently smooth 
that the (Aq,A\) phase transition is controlled by F(gk,gk), c.f. (4.10) and (4.26), and is therefore 
3 rd order. It is plausible that these results, valid in any finite order of the double perturbation 
theory, should also hold at a strong coupling A satisfying Nf/N < 1/A as well. Gravitational 
approximations are valid at large g s N, so that the backreaction on D branes is small. To keep 
open string perturbation theory in control, one also needs g s Nf <C 1 which means Nf/N <C 1/A in 
the large N limit. We conjecture that the third order phase transition is valid if these conditions 
are meet. 2 This would then imply that the string theory dual to the gauge theory would exhibit 
a 3 rd order Hawking-Page transition (dual to the confinement /deconfinement transition) in the 
supergravity limit of the string theory. No such transition has yet been found when (Nf/N) = 
finite in the large N limit, as one must go beyond the probe approximation, and consider the 

2 We thank the referee for this argument. 



9i = 91 - 1/2 • 



(4.25) 



S eff = - F (9k,9k) + f(9k,9k) 



(4.26) 
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back-reaction of the flavor branes. There have been preliminary studies involving the back-reaction 
of flavor branes [10, 11], but are not able to discuss the Hawking-Page phase transition from a 
thermal ensemble to a large black hole. Discovery of such a 3 rd order transition in the string dual 
to a large N gauge theory with (Nf/N) / 0, would be an exciting confirmation of duality in the 
context of the ideas of this paper. 
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